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Many neurons exhibit membrane potential resonance 
(MPR), a peak in the membrane impedance amplitude 
(|Z|) in response to oscillatory inputs at nonzero fre- 
quency (fmax) [!]• MPR arises from nonlinearity and 
timescales of voltage-gated currents and may set fre- 
quency of network oscillations. Pacemaker PD neurons 
of the crab pyloric network show MPR whose fmax is 
correlated with the network frequency (~ IHz) [2]. In 
contrast, the LP follower neuron shows a higher /^^x of 
~ 1.4 Hz. The impedance profile of biological PD and 
LP neurons and the model neuron was measured using 
a logarithmic ZAP function (/miH=0.1 Hz, /„ax=4 Hz) in 
voltage clamp (V;o„,=-60mV and V/,/^/,=-30mV). The 
/msA^in biological PD neurons increases if either V/o^ or 
Vhigh are increased [3], whereas the LP newon f^ax is 
only sensitive to V/j/^/,. Additionally MPR in the PD neu- 
rons is sensitive to blockers of Ica and //,. We hypothe- 
size that: (1) many combinations of parameters can 
produce MPR in PD and LP neurons; (2) The MPR 
mechanism in LP is distinct from PD. 

Experimentally, /cais difficult to measureand therefore 
a top-down approach is adopted to elucidate the contri- 
butions of /caand 1^ to MPR in PD and LP. Because 
resonance depends on the kinetics of /caand Ih, a brute- 
force sampling of the parameter space is computation- 
ally unfeasible and, therefore, we search for model para- 
meters using a genetic algorithm. The biological data 
were used to constrain the range of leak, Ica and 1^ 
parameters in a single-compartment model. The genetic 
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algorithm, NSGA-II [4] was used to optimize the MPR 
profile and produce a population of optimal models. A 
sensitivity analysis of MPR attributes on model para- 
meters was done in these models. 
The distributions of optimal parameters were tightly 

constrained for gi^^i,, V'y, Ca_act. ^M._Ca_inact and T Cajnact- 

Additionally, strong correlations were observed between 
T_ca_act and T ca_inact (negative), between Vy, ca_act and 
V^/, Cajnact and between ^c^and VV, ca_act (negative). In 
models with low Ihtfmax correlated strongly with the fre- 
quency which /(ja peaked, which is controlled by x a^t 
and T Cajnact- The parameter sensitivities also support the 
sensitivity to /catime constants, demonstrating potential 
targets for neuromodulation. 

The MOEA was also used to optimize the /,„axshifts 
with V/o^ and V},ig/, to produce two model groups with 
properties that correspond to the differences between 
PD and LP. These results suggest that ^a;tshift is due to 
different activation rates of Ih and therefore these two 
neurons may generate MPR through different mechan- 
isms; a result which we aim to test experimentally. 

Many neurons display emergent properties in response 
to oscillatory inputs, such as amplified responses in cer- 
tain frequency bands. These properties may be important 
in shaping coherent network activity. The underlying 
nonlinearities and time scales that shape specific features 
of impedance profiles can be used to link sub-threshold 
dynamics to supra-threshold voltage responses. We have 
used an MOEA to understand the multiple underlying 
ionic mechanisms that generate resonance and explained 
how PD, and not LP,/,„aa:Can be adjusted according to 
different input amplitudes. 
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